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Abstract: We study real time correlators in strongly coupled M = 4 supersym- 
metric Yang-Mills theory on dS$ x S 1 , with antiperiodic boundary conditions for 
fermions on the circle. When the circle radius is larger than a critical value, the dual 
geometry is the so-called "topological AdS$ black hole" . Applying the Son-Starinets 
recipe in this background we compute retarded glueball propagators which exhibit 
an infinite set of poles yielding the quasinormal frequencies of the topological black 
hole. The imaginary parts of the propagators exhibit thermal effects associated with 
the Gibbons-Hawking temperature due to the cosmological horizon of the de Sitter 
boundary. We also obtain R-current correlators and find that after accounting for a 
small subtlety, the Son-Starinets prescription yields the retarded Green's functions. 
The correlators do not display diffusive behaviour at late times. Below the critical 
value of the circle radius, the topological black hole decays to the AdS$ "bubble of 
nothing" . Using a high frequency WKB approximation, we show that glueball cor- 
relators in this phase exhibit poles on the real axis. The tunnelling from the black 
hole to the bubble is interpreted as a hadronization transition. 



1. Introduction 



Time dependent backgrounds in gravity and in string theory are of great interest 
from the standpoint of the AdS/CFT correspondence [p], § and related holographic 
dualities between gauge theories and gravity. Time dependent classical gravity 
backgrounds, in asymptotically Anti-de-Sitter spacetimes, can potentially provide 
a fully nonperturbative description of non-equilibrium phenomena in the strongly 
coupled dual gauge theories. Such non-equilibrium physics in field theories arises, 
most notably, in cosmology and in heavy ion collisions at RHIC. To understand how 
gauge/gravity dualities work for such processes, it is important to investigate how 
holography applies in various examples with explicit time dependence. In this paper 
we attempt the holographic computation of real time correlators of the boundary 
gauge theory dual to the time dependent, asymptotically locally AdS backgrounds 
found in |, 

The authors of || studied the double analytic continuations of vacuum solu- 
tions such as Schwarzschild and Kerr spacetimes providing examples of smooth, time 
dependent solutions called "bubbles of nothing" || [10], [11]. These asymptotically 
flat solutions were generalized to asymptotically locally AdS spacetimes in |3], |], 
by considering the double analytic continuations of AdS black holes 1 . The bubbles 
are obtained by analytically continuing the time coordinate to Euclidean signature 
t — > %x where x is periodically identified, and a polar angle 9 — > it. In addition, the 
X circle has supersymmetry breaking boundary conditions for fermions. The result- 
ing "bubbles" undergo exponential de Sitter expansion (and contraction). For the 
asymptotically locally AdS$ x S 5 case the conformal boundary of the geometry 
is dSs x S 1 . The corresponding dual field theory, Af = 4 SYM, is thus formulated on 
dSa x S 1 with antiperiodic boundary conditions for the fermions around S 1 . Each 
of the two AdS-Schwarzschild black holes (the small and big black holes) yield an 
AdS bubble of nothing solution, only one of which is stable. The bubble of nothing 
geometries are vacuum solutions with cosmological horizons M and particle creation 
effects. 



It was realized in |5|, |6|, |7J that there is another spacetime with the same AdS 
asymptotics as the bubble geometries, with dS^ x S 1 conformal boundary. This is 
the so-called "topological black hole" 2 - a quotient of AdS space obtained by an 
identification of global AdS§ along a boost [17, 18|. It is the five dimensional analog 

1 For the classifications of solutions obtained by analytically continuing black hole solutions, see 

El- 

2 The term "topological AdS black hole" has also been used to refer to black holes with a hyper- 
bolic horizon having a non-trivial topology. In the AdS/CFT context these have been studied in 
H , |l5| , [lq] and references therein. 
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of the BTZ black hole [19, 2C]. The topological AdS black hole can also be obtained 
by a Wick rotation of thermal AdS space. As Euclidean thermal AdS space can 
be unstable to decay to the big AdS black hole via the first order Hawking-Page 



transition |2T], 22], a similar instability is associated to the topological AdS black 
hole. In this case the topological AdS black hole is unstable to semiclassical decay 
via the nucleation of an AdS bubble of nothing. The associated bounce solution is the 
Euclidean small AdS-Schwarzschild black hole which has a non-conformal negative 
mode. The topological black hole becomes unstable only when the radius of the 
spatial circle becomes smaller than a critical value (in the Euclidean thermal setup 
this is when the temperature exceeds a critical value). Precisely such an instability 
to decay to "nothing" was, of course, first noted for flat space times a circle having 
antiperiodic boundary conditions for fermions M. 



The two different geometries described above are dual to two different phases 
of strongly coupled, large N gauge theory formulated on dS 3 x S 1 . As in the usual 
thermal interpretation wherein the field theory lives on S 3 x S 1 , the two phases 
are distinguished by the expectation value of the Wilson loop around the S 1 . In 
the bubble of nothing phase, the circle shrinks to zero size in the interior of the 
geometry and the Wilson loop is non-zero, indicating the spontaneous breaking of 
the Z/v symmetry of the gauge theory. The topological black hole phase is Z^v 
invariant. Unlike the thermal situation however, the spontaneous breaking of Zjv 
invariance is not a deconfinement transition since the circle is a spatial direction and 
not the thermal circle. 



Our primary motivation in this article is to understand how the behaviour of real 
time correlators in the two geometries reflects the properties and distinguishes the 
two phases of the M = 4 theory on dS 3 x S 1 . Since the de Sitter boundary has its own 
cosmological horizon accompanied by a Gibbons-Hawking radiation |23| , this should 
also be reflected in the properties of the boundary correlation functions. An interest- 
ing feature of both the geometries in question is that infinity is connected, i.e. the 
asymptotics is unlike the (AdS-)Schwarzschild black hole whose asymptotics consists 
of two disconnected boundaries. This means the Schwinger-Keldysh approach in p4 



is not directly applicable. It would be interesting to understand how to apply that 
idea and also the recently proposed prescription of J25L BB1 in the present context. 



Instead we simply use the Son-Starinets prescription |27|, ^] to compute real 
time correlators in the topological AdS black hole geometry, by requiring infalling 
boundary conditions at the horizon of the black hole. 



We find that in the topological black hole phase, retarded scalar glueball correla- 
tors (homogeneous on spatial S 2 slices of dS 3 ) have a simple description in frequency 
space. They have an infinite number of poles in the lower half of the complex fre- 
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quency plane. As in the case of the BTZ black hole and other well known examples, 
these poles represent the black hole quasinormal frequencies [p7 ]. The Green's func- 
tions have imaginary parts and display features closely resembling thermal physics. 
These features are naturally associated to the Gibbons-Hawking temperature due to 
the cosmological horizon of de Sitter space. This suggests that the Af = 4 theory 
on dSz x S 1 is in a plasma-like or deconfined state in the exponentially expanding 
universe. 



We further investigate real time correlators involving spatial spherical harmonics 
of conserved R-currents to find whether they exhibit transport properties, i.e., if 
they relax via diffusion on the expanding spatial S 2 slices of dS^. Applying the 
Son-Starinets recipe (here we have to acount for a certain subtlety involving discrete 
normalizable mode functions in de Sitter space) we find that the retarded propagator 
of the R-current does not appear to relax hydrodynamiccally. This is likely due to the 
"rapid" expansion of de Sitter space, the expansion rate of dS^ being of the same order 
as the Gibbons-Hawking temperature. The real time correlators are represented in 
the form of a de Sitter mode expansion, which allows to identify a natural frequency 
space correlator. This latter object has isolated poles in the lower half plane and at 
the origin, and its imaginary part exhibits the features characteristic of a thermal 
state. 



When the spatial circle is small (relative to the radius of curvature of dS^), below 
a critical value, the topological black hole decays into the AdS bubble of nothing. In 
this geometry, correlation functions are not analytically calculable. However, scalar 
glueball propagators can be calculated in a WKB approximation. We show that in 
this approximation, the correlation functions have an infinite set of isolated poles 
on the real axis in the frequency plane. We interpret this naturally as high mass 
glueball-like bound states of the field theory. The transition from the topological 
black hole to the bubble of nothing by tunelling is interpreted as a hadronization 
process. A related picture of hadronization was discussed in |HJ . 

The paper is organized as follows. In Section 2, we review properties of the 
topological AdS black hole. In Section 3, we perform the detailed holographic com- 
putation of retarded propagators of the spatially homogeneous scalar glueball fields. 
We also calculate R-current correlation functions. Section 4 is devoted to a WKB 
analysis of Green's functions in the bubble of nothing phase. We summarize our 
results in Section 5. 
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2. The Topological AdS Black Hole 



The so-called topological black hole of 0, [18| in AdS 5 is an orbifold of AdS space, 



obtained by an identification of points along the orbit of a Killing vector 

Z = -p^(x A d 5 + x 5 d A ), (2.1) 

-KAdS 

where r x is an arbitrary real number and the AdS space is described as the universal 
covering of the hyper-surface 

— x Q + x 1 + x 2 + x 3 + a; 4 — x 5 = — R^si (2-2) 

-^Ads being the AdS radius. In Kruskal-like coordinates which cover the whole space- 
time, the metric has the form 

d °' = <WV + §3*$ r x* 2 (2.3) 

where x is a periodic coordinate with period 2ir. The four coordinates y fl , (fi = 
0, . . .3) are non-compact with the Lorentzian norm y 2 = y >1 y u i]^ L/ such that — 1 < 
y 2 < 1. Locally, the spacetime is anti-de Sitter with a periodic identification of the 
X coordinate, 

X~X + 2vr. (2.4) 

The conformal boundary of the spacetime is approached as y 2 — > 1, and it is dS^ x S . 
The boundary conformal field theory is therefore formulated on a three dimensional 
de Sitter space with radius of curvature i?AdS times a spatial circle of radius r x . 

The geometry has a horizon at y 2 = 0, which is the three dimensional hypercone, 

vl = yl + yl + yl (2.5) 

and a singularity at y 2 = —1. The singularity appears because the region where 
the Killing vector has negative norm needs to be excised from the physical space- 
time to eliminate closed timelike curves. The hyperboloid y 2 = — 1 is a singularity 
since timelike geodesies end there and the Killing vector d x generating the orbifold 
identification has vanishing norm at y 2 = —1. The topology of the spacetime is 
M 3 ' 1 x S 1 , in contrast to that of the AdS-Schwarzschild black hole which has the 
topology R 1 ' 1 x S 3 . For this reason, infinity is connected in this geometry unlike in 
the usual Schwarzschild black hole which has two disconnected asymptotic regions. 

Finally, it is possible to rewrite the metric in Schwarzschild-like coordinates by 
introducing 

y2 =i>- - =coth f^)' TrlL (2-6) 

^ y V^Ads/ #AdS (l + Vo-Y 2 ) 2 



i=l 
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Figure 1: The global structure of the topological AdS black hole spacetime. The singularity is 
the hyperboloid — j/q + UiVi = ~~ 1 and the horizon is at the cone j/q = ViVi- 

These coordinates only cover the region y 2 > which is the exterior of the topological 
black hole. Locally, the metric takes the form (simply related to Eq.(ll) of after 
a coordinate transformation), 

ds 2 = 

^vfks + (^) 2{r2+ Rl ^ 2 + r2 (-£ + cosh2 {£) ^ 

(2.7) 

The Euclidean continuation of the metric yields thermal AdS space due to periodicity 
of the x coordinate. Hence, the topological black hole metric (exterior to the horizon) 
can also be obtained following a double Wick rotation of global AdS spacetime and a 
periodic identification of the x coordinate. In the Schwarzschild-like coordinates, the 
horizon of the topological black hole is at r = 0. It is clear that each slice of constant 
r is a dSs x S 1 geometry. This metric, while locally describing AdS space, differs 
from it globally due to the identification x ~ X + 2vr. Note also that the spatial S 1 
remains finite sized at the horizon, with radius Rs^ = r x . 

It is interesting to see that we can get a better understanding of the geometry 
in the vicinity of the singularity at y 2 = — 1 by zooming in on the the metric ( |2.3| ) 
in this region. Introducing the coordinates 

yo = (1 -5) coshe, Y = (1 -5) smhe < 5 < 1, (2.8) 
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we find 

ds 2 ss R 2 AdS {-d5 2 + de 2 + sinh 2 e dQ 2 2 ) + r 2 J 2 d X 2 , (2.9) 
which is the metric for Milne spacetime (in the 5, x directions). 

3. Real time correlators in the Topological AdS Black Hole 

We will compute real time correlators in the Yang-Mills theory on the boundary of 
the topological AdS 5 black hole (TBH) following the recipe of Son and Starinets [f27| 
in the Schwarzschild-like patch ( |2.7| ) of the black hole. Viewing the topological black 
hole as a Wick rotation of (Euclidean) thermal AdS space, one expects that such 
correlators can also be obtained by an appropriate analytic continuation of Euclidean 
Yang- Mills correlators on S 3 x S 1 in the confined phase (the Z^r symmetric phase) 
with anti-periodic boundary conditions for fermions. Since the relevant Wick rotation 
turns the polar angle on S 3 into the time coordinate of de Sitter space, a complete 
knowledge of the angular dependence of Euclidean correlators on S 3 x S 11 would be 
necessary. However, finite temperature Yang-Mills correlators on S 3 and at strong 
coupling, have not been calculated explicitly, so we will not follow the route of 
analytic continuation. Instead we will directly calculate the real time correlators 
using the holographic prescription of Son and Starinets applied to the topological 
AdS black hole geometry. 



3.1 Scalar wave equation in the Topological Black Hole 



To extract field theory correlators, we first need to look for solutions to the wave 
equation in the region exterior to the horizon of the topological black hole. It is 
instructive to write the metric for the black hole in the Schwarzschild form of [|7| 



ds 2 = R 



AdS 



dp' 



(P 2 - 



i) + G&) 2 p2dx2 + {p2 _ i} ( ~ dr2 + cosh2 r dnl) 



(3.1) 



where we have introduced the dimensionless variables 

P = 



r/R 



AdS) 



+ 1, 



T 



R\dS 



(3.2) 



The conformal boundary of the space is approached as p —>■ oo while the horizon 
is at p = 1, where the coefficient of dr 2 vanishes. The slices with constant p are 
manifestly dS^ x S 1 . 
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The scalar fields in this geometry have a natural expansion in terms of harmonics 
on the S 2 x S 1 spatial slices 



/d 
. 



>,r). (3.3) 



£,m,n 

The normal mode expansion above involves spherical harmonics on S 2 , the discrete 
Fourier modes on S 1 and %(u, r) which solve the scalar wave equation on dS^. 

—^d T (cosh 2 r d T T t {u, r)) + £ -^±^-T e (u, r) = -[v 2 + l)T t (u, r). (3.4) 
cosh r cosh T 

For every i e Z, the equation has two kinds of solutions that will be relevant for 
us: i) normalizable modes labelled by integers —iu = 1, 2, . . . £; and ii) delta-function 
normalizable modes labelled by a continuous frequency variable u 6 R. We will 
return to this point when we discuss R-current correlators. General solutions to this 
equation can be expressed in terms of associated Legendre functions 

T t (u, r) = ( Ai Pf(tanhr) + B x Qf (tanhr)) . (3.5) 

coshr 

In the usual approach to quantizing free scalar fields in de Sitter space, the integration 
constants Ai and Bi are determined by the choice of de Sitter vacuum |3~I|, [321 



33|| . However, in the present context, the constants will be specified by picking 
out inf ailing wave solutions at the horizon of the topological black hole. These are 
the holographic boundary conditions relevant for real time response functions in the 
strongly coupled field theory on dS 3 x S 1 . 

It is useful to see the scalar wave equation in this background recast as a 
Schrodinger equation, using Regge- Wheeler type variables 

M = -ln(^-i — ] or p = cothw, (3.6) 



2 VP" 1 , 
and 

®n = VP(P 2 ~ 1) $n, (3.7) 

$ being the scalar field in the bulk. In these coordinates, the horizon is approached 
as it — > oo while the conformal boundary is at u = 0. Following the above coordinate 
and field redefinitions, the Schrodinger wave equation in the TBH geometry is 

d 2 

- -j-p ^n(y, u) + V n (u) V n (v, u) = v 2 V n (v, u), (3.8a) 

V n (u) = ((MR AdS ) 2 + f )^-L- + (| + %^n 2 )— L- . (3.8b) 

smh u * cosh u 
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Figure 2: The Schrodinger potential for the topological AdS black hole. 



Here, \l/ n = $ n \/ p(p 2 — 1) and we have allowed for a generic non-zero mass M since 
we will eventually be interested both in the massless and massive cases. As expected 
for AdS black holes, the potential decays exponentially near the horizon u — > oo, 
while blowing up near the boundary at u — > 0. In the near horizon region where 
the potential vanishes, the solutions with v > are travelling waves and there is a 
natural choice of incoming and outgoing plane wave solutions. For any n and M, the 
equation has analytically tractable solutions in terms of hypergeometric functions. 
We will use these to calculate the retarded Green's functions for the boundary gauge 
theory (i.e., Af — 4 SYM) at strong coupling on dS 3 x S 1 . 

Although analytical solutions exist for all non-zero n and M, we will restrict 
attention, for simplicity, to two special cases: i)n^0 and MR^ds = 0; ii) n = and 
MR^ds 7^ 0. In each of these two cases the radial equation is solved by two linearly 
independent hypergeometric functions: 



M = 0; n^O : 

( p 2 _ ^-if-i 2Fi _ + n Ems^ §-i(i/ + n ); 1 - m p 2 ) 

+ C 2 p^AdsAx x 

(p 2 - 2 F X (-l-i(v- n ), § - - n ); 1 + m ; p 2 ) . 

(3.9) 
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Similarly, for massive fields with n — 0, the two linearly independent solutions are 
M/0; n = : 

$o(f,p) = 

d (p 2 - l)-^ 4 -^ 2 Fx (|(3 - A) - |(3 - A) + 3 - A; -(p 2 - l)" 1 ) 
+ C 2 (p 2 - l)-i A 2 F X (|(A - 1) - ii„, |(A - 1) + ife; A - 1; -(p 2 - l)- 1 ) 



A = 2 + v/4+(Mi? AdS ) 2 . 

(3.10) 

The correct linear combination, relevant for the holographic computation of correla- 
tors, is picked by applying the requirement of purely infalling waves at the horizon 
of the topological black hole at p ~ 1. 



3.2 Scalar glueball correlator 



As argued in @, [7|, the topological black hole in AdS space is automatically a so- 
lution to the Type IIB supergravity equations of motion, since it can be obtained 
via a double Wick rotation (and an identification) of AdS§ x S 5 . The M = 4 super- 
symmetric Yang-Mills theory on the dS^ x S 1 conformal boundary of the topological 
black hole, has two 5*0(6) singlet, scalar glueball fields 

Q{x, t) = TrF^F^; G(x, t) = TrF^F^. (3.11) 



These are dual to the dilaton and the RR-scalar in the Type IIB theory on the bulk 
spacetime and both solve the massless (corresponding to A = 4 operators on the 
boundary) scalar wave equation in the topological black hole geometry. The retarded 
propagators for the scalar glueball fields are known on IR 3,1 at weak coupling both at 
zero and finite temperature [[53J]. Since the operator is chiral primary in the Af = 4 
theory, at zero temperature its propagator on M 3,1 receives no quantum corrections 
and the strong coupling results from supergravity are in exact agreement with those 
of the free field theory. At finite temperature, however, when supersymmetry is 
broken, strong and weak coupling results on M 3 ' 1 differ [^7, 34]. Computations of the 
glueball correlator also exist in the free M = 4 theory at finite temperature and on 
a spatial S* 3 , both in the confined and deconfined phases |34|]. Their strong coupling 
counterparts have not been determined. 
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The present case, with the field theory on dS$ x S 1 , is intriguing for the following 
reasons. First, there is the lack of supersymmetry, due to antiperiodic boundary 
conditions for fermions around the spatial S 1 . Secondly, the boundary field theory 
sees a cosmological horizon on dS^ accompanied by its thermal bath. It would be 
interesting to observe the emergence of the boundary Gibbons-Hawking temperature 
from a holographic calculation of its correlators at strong coupling. Finally, when 
the radius of the boundary S 1 decreases below a critical value, r x /i?Ads < 5775; the 
topological black hole decays via a bounce to the small "AdS bubble of nothing" . We 
would like to understand how boundary field theory correlators at strong coupling on 
dSs x S l change across this transition. The transition from the topological black hole 
to the Bubble of Nothing is a breaking transition. This is understood precisely as 
in the Euclidean (finite temperature) situation, due to a non-zero expectation value 
for the Wilson loop around the spatial S 1 . 

Curiously, it is apparent that in the classical supergravity approximation, the 
bulk scalar glueball correlators are insensitive to fermions and their boundary con- 
ditions around the spatial S 1 . It would be interesting to understand whether this is 



related to large- iV volume independence J35], [36| in the symmetric phase. 



We are primarily interested in real time response and for the sake of simplicity, 
we will first study only the response functions for glueball fluctuations that are 
homogeneous on the spatial S 2 slices at the boundary, i.e., 

G r (t,t' ;n;l = 0) = 

dndtf fd X /wr „ /0 N „ /0/n ... (3-12) 



— 1 



/SIT / |^""e(r-r')([6(!),x,r), S(ff,0,x<) ] ). 



We will work with the dimensionless variables r = t/Rj^, t' = t'/RAdS an d restore 
appropriate dimensions when necessary. As we will see when we look at R-current 
correlators, it is straightforward to generalize to inhomogeoneous fluctuations on the 
spatial sphere. 

For the moment we focus attention on the s-wave (£ = 0) retarded correlation 
function of the scalar glueball operator. Also, for the s-waves, the correlator turns 
out to be a function of (r — r') so that it is natural to define the temporal Fourier 
transform of this as, 

/oo 
dr e - iu{r - T,) G r (t, T';n;£ = 0). (3.13) 
-00 

To calculate it at strong coupling and in the large radius regime (r x > -f^f), we solve 
the dilaton wave equation which is the equation for a massless, minimally coupled 
scalar field in the background of the topological AdS black hole. 
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3.2.1 Spatially homogeneous case with n = £ = 



We begin by looking at the spatially homogeneous response functions, I = n = 0, on 
the dSs x S 1 slices. The solutions to the radial part of the Klein-Gordon equation in 
the massless limit are the hypergeometric functions 



and 



*< 2) (M = (P 2 - l)' 2 2F1 ^ , ^ ; 3 ; . (3.15) 



For the Z = modes, the temporal dependence is also particularly simple, and has a 
natural interpretation in terms of positive and negative frequency states 

Solving the Dirichlet problem and extracting correlation functions holographically 
requires us to first pick the correct linear combination of the two solutions which 
is smooth near the horizon p — > 1 and represents an incoming wave falling into the 
horizon. In the near horizon region, the asymptotic form of the solutions is: 

, i-* „ r(-«/)r(^) Wv r(iv)r(2=p) 

i (2(p - eT" l r( ^ } 2 ^ + i (2( P - l))"" *" L ^2 ; 

(3.17) 

and 

•FW-l)- (2(p-l))-— (2^-1))-—^^,. (3.18) 



Note that these modes diverge as l/y/(p — 1) near the horizon. However, em- 
ploying the measure implied by the bulk metric y^g ~ (p 2 — l) 5 ^ 2 , these are still 
normalizable in the vicinity of the horizon. We can now pick a linear combination 
such that only incoming positive frequency waves are allowed at the horizon. This 
means, assuming that are the positive frequency modes with Re(z/) > 0, the 
solution to the radial wave equation should behave like (p— 1) _ near the horizon. 
In conjunction with this we have the properly normalized boundary behaviour as 
P -> 00, 

V, P - 00) - 1 + . . . ; <4 2 V, p - 00) - ^ + • • • • (3-19) 
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The complete solution to the boundary value problem for a massless scalar with 
I = n = 0, in the topological AdS black hole is then, 



(3.20) 



Following the holographic prescription |27, 28| for computing real time correlators, 
the Yang-Mills retarded correlation function is obtained by analyzing the boundary 
terms from the on-shell scalar action 
N 2 

- I (It I <n> I <l\ ( ,'''\, -a 4>( r. o) DAM r. o) 

I p— >oo 



16tt- 



s 



(3.21) 



where, for the spatial s-wave we have defined 

/OO i 
-oo 2vr 



(3.22) 



Putting together the explicit expressions for T^iv^ r), and the boundary behaviour of 
the solution fl3.2U| ) we are immediately led to the (unrenormalized) s-wave retarded 
correlator in frequency space, including all contact terms (finite polynomials in the 
frequency v) 



Gr(u;0) 



N 2 



"1(1 + - 2 ) 2 



3 — izv 



if, 



3 + iv 



in coth (n\(v + i)) 



+ 1(1 + v 2 ) 2 (Inp - 1E + 1) \ p _ + 1(1 + u 2 )p 2 \ p _) 



(3.23) 



We remark that unlike the case of the Poincare' patch description of AdS space, the 
non-normalizable solution in the topological AdS black hole (akin to global AdS), 
contains a term proportional to 1/p 2 in its near-boundary asymptotics, but it only 
contributes a quadratically divergent contact term in the correlation function above. 



The divergent and scheme-dependent contact terms can be minimally subtracted 
away to yield the renormalized, retarded Green's function. Up to now we have been 
working with dimensionless variables, corresponding to a de Sitter boundary of unit 
curvature. Restoring dimensionful constants with the replacement 

v -> vR, (3.24) 

R being the radius of curvature of the boundary dS$ (or the inverse Hubble constant), 
the renormalized retarded Green's function continued into the complex frequency 
plane is 

N 2 



3 — ivR 



HvR 

;i + v 2 r 2 ) 



(3.25) 
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Figure 3: The analytic structure in complex frequency plane, of the massless, spatially homo- 
geneous (I = n = 0), retarded Green's function in the Zjv symmetric phase, corresponding to the 
topological AdS black hole. 



For Im(z/) = 0, its real and imaginary parts match ( 3.23| ), and the function is analytic 
in the upper half plane, with only isolated simple poles in the lower half plane at 



u k = -i(3 + 2k) 



R 



(3.26) 



As argued in |27[ , poles of the retarded correlator in a black hole background coincide 
with the quasinormal frequencies of the black hole. The quasinormal frequencies and 
the retarded glueball correlator found above, for the topological black hole in AdS 5 , 



are remarkably similar to the corresponding objects in the BTZ black hole [27 . 



3.2.2 Non-zero momentum along S 1 and I = 



It is relatively easy to allow for a non-zero discrete momentum n/r x along the spatial 
S 1 . This requires the modes ( |3.9| ), to solve the boundary value problem in the topo- 
logical black hole background. Following the same steps as in the s-wave correlator, 
(after tedious algebra) we find that the retarded Green's function (with dimensionful 
constants restored) is 



N 2 



1287T 5 



2*fl("+£) 



V 



^) 2 + R~ 



-2 



X 



neZ. 



2iR(v - 

r x 

/ - ^) 2 R 2 + 1 



(3.27) 



13 



When n — 0, this matches our expression for the s-wave correlator (|3.25| ). The 
Green's function has nonanalyticities only in the lower half plane, with simple poles 

at 

v± = -i(3 + 2k)R~ 1 ±— ; k,neZ, (3.28) 



giving the quasinormal frequencies of the topological black hole, with non-zero mo- 
mentum along the spatial S 1 . Interestingly each simple pole at n = 'splits' into two 
simple poles at non-zero n. Our expressions for the correlation functions on dS% x S* 1 



n/r^ 

- i 3/R • 
-i7/R • 
4 5/R • 
O 
O 



Figure 4: Simple poles in the frequency plane, of the massless retarded Green's function with 
non-zero momentum along the spatial circle. 

in the Z^r symmetric phase satisfy the obvious consistency check - in the high fre- 
quency/large momentum limit uR,n 3> 1, they reproduce the flat space Green's 
function for the scalar glueball operator 

_ N 2 

i/R,n>l 1287T 2 



G R (u;n) 



W - 



n 2 \2 



ln(z/ 2 - S 



(3.29) 



3.3 Thermal effects and the Gibbons-Hawking temperature 



De Sitter space has a cosmological horizon and an associated Gibbons-Hawking tem- 
perature PB| . We therefore expect our boundary (dS 3 x S 1 ) field theory correlators 
to exhibit thermal properties. For real frequencies z/, the digamma functions have 
an imaginary part so that 

N 2 



Im GrU 1 ; n 



Im(i/)=0 2567T 



v — 



n \2 



+ R- 



y+^f + R- 2 ) x 



ttR 

coth( — 



%R 



TiR 



+ coth — (u+^-iR 



-v 



(3.30) 
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Here we have used tanh(:r) = coth(x + i|), to cast the result in a form that will make 
the connection to thermal physics explicit. 



In flat space and in free field theory at finite temperature T / 0, the scalar 
glueball propagator , with zero spatial momentum and frequency uj is proportional 
to the digamma function p| 



^Mriatspace = "^2 u ^ (4^) + analytic. (3.31) 
The imaginary part of the flat space glueball correlator is the spectral function 

Im G«M| Flatspace = -^Vooth (^) . (3.32) 

In weakly coupled, or free field theories at finite temperature (on flat space), the one- 
loop spectral function reflects the physical effect of 'Bose enhancement', following 
from stimulated emission of bosons into the heat bath. In perturbative field theory 
on flat space, this manifests itself as an enhancement of the decay rate of an unstable 
boson in a heat bath, at rest with energy u, by a factor (relative to the vacuum decay 
rate) 

coth (ir) = 1 + 2nB ® = 1 + ( 3 - 33 ) 

While there may not be an obvious way define a spectral representation in de Sitter 
space, the similarity between our strongly coupled de Sitter space result ( |3.30| ) and 
( |3.32| ) is obvious. In particular, it allows the identification of a temperature in de 
Sitter space 

r„ = ^ (3.34) 

which is precisely the value of the Gibbons-Hawking temperature. Note that despite 
the similarity between the expressions for dS% and thermal correlators in flat space, 
there is a crucial difference - the frequency or 'energy' appearing in the Bose-Einstein- 
like distribution function in de Sitter space, is not the real frequency v ( |3.30| ), but 



in fact v — iR . This difference can be traced back to the definition of our positive 
and negative frequency modes ( |3.51 ). For real v, the positive frequency modes are 



red-shifted away in the far future. To get propagating modes in the future, we would 
need to choose v = uo + iR^ 1 with w6l 



It is interesting to note that our results for the retarded glueball propagator in 
the strongly coupled field theory (in the Zn symmetric phase) on de Sitter space 
closely match one- loop weakly coupled field theory calculations |37[| . It would be a 
straightforward calculation to check whether there is exact agreement between weak 
and strong coupling results on dS 3 x S 1 . We leave this excercise for the future. For 
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now, we only make the following observation, which suggests that the scalar glueball 
correlator in the Zat symmetric phase should not be renormalized. 



It has been argued in |35|, |36| that correlation functions of large N gauge theories 
in the TL^ symmetric phase, with some or all spacetime directions compactified, 
should be independent of the volume of the compact directions. In the present 
situation this would imply that on dSz x S 1 with antiperiodic boundary conditions for 
the fermions, large N correlators in the Z N symmetric phase (r x > RAds/^V^) should 
be independent of the radius of the S 1 . In particular then, for perturbations which 
are homogeneous along the circle, the correlation functions should be independent 
of r x and should match up with the result on dS$ x R. The latter is obtained by 
a (double) Wick rotation of S 3 x R. Since TrF 2 is a chiral primary in the Af = 4 
theory and its correlator on S 3 x R is not renormalized by interactions, one would 
expect this to be true also on dS^ x R. 



3.4 The massive case 



The holographic calculation of correlation functions in the topological AdS black hole 
can be easily extended to massive scalars. In the context of the Type IIB theory, 
such massive states are stringy excitations with masses 

scalar field of mass M in the bulk is dual to a scalar operator Oa in the field theory 
with conformal dimension A = 2 + a/4 + (Mi?Ads) 2 - We will study below the free 
massive scalar in the bulk geometry to extract information on the analytic structure 
of correlators of high dimension operators in the field theory. 



There are two primary motivations for looking at high dimension operators in the 
field theory: i) The works of |B8L BD| have demonstrated that propagators of heavy 



fields, in the geodesic approximation, may be used to probe the bulk geometry behind 
horizons and perhaps extract information on singularities behind such horizons, ii) 
One of our main goals is to look for signatures of the transition between a Z^r 
symmetric phase and a Z^r broken phase. In the bulk theory, the latter phase is the 
small bubble-of-nothing geometry. Correlators in this latter geometry can only be 
computed using an eikonal (WKB) approximation involving high frequencies and/or 
large masses. 



Extending the holographic analysis done above for massless fields, to massive 
scalars in the topological black hole geometry, we find that the frequency space 
correlator is 

^)^ A r ^ (A - 1 -^ ) ) 2 2 r(3 - A) , (3.35) 

r(i(3-A- Wj R)) 2 r(A-i)' 
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where the normalization constant Ca = 2(A — 2)e 2< > A A \ with e — >• as the boundary 
is approached. 



Double poles 



Figure 5: Double poles in the lower half plane for the massive (retarded) propagator at zero 
spatial momentum. 



In the massless limit MR — > 0, this reproduces the expression (|3.25|) found 



previously, after subtracting an additional divergent contact term. The massive 
correlator has an analytic structure that is qualitatively distinct from the massless 
case. In particular the retarded correlator has an infinite set of double poles and 
simple poles at 

u k = -i(A-l + 2Jfe)4; keZ. (3.36) 
R 

The significance of the appearance of double poles in the massive retarded propagator 
is not entirely clear. Such double poles have also appeared in 2d CFT correlators 
with non-integer conformal dimensions from the BTZ black hole |]27|] , at zero spatial 
momentum. At finite spatial momentum (e.g. n ^ 0) we expect the double poles to 
split into simple poles. 

For real frequencies, the massive correlator also has an imaginary part 

Tm r <v\ - *V r(3 ~ A) sin(7rA) sinh(TrzAR) 

Im G R {v) - -_£ A __ | r( l (3 _ A _^ ))cos(f(A _^ ))r ( 3 - 37 ) 

The (de Sitter) thermal origin of this result is not as explicit as for the massless scalar. 
However, after identifying the de Sitter Hawking temperature to be T H = R~ 1 /2tt, it 
is worth comparing the above expression with the imaginary part of the propagator 
in two dimensions for large non-integer conformal dimension deduced from the non- 



extremal BTZ black hole p7 |. The similarities between the two results, particularly 



the numerator of ( |3.37| ), are striking. 
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The large mass, high frequency limit of this result is easily obtained, using Stir- 
ling's approximation 

(3.38) 



'2tt^= e~ z z z . 

<z 



When the masses are taken to be large so that MR ^> 1, then A « MR. In this 
high frequency, large mass limit it is useful to define a rescaled frequency variable 



v = 



v 



vR, MR > 1, 



so that 



G R {v) ~ C A 



■ ~\ MR(1- 
IV x 



■1 — IV 



MR(l+iv) 



(3.39) 



(3.40) 



Here we have ignored an overall (real) phase due to frequency independent coefficients 
in the large mass limit. At first sight, a potentially problematic feature of this 
approximation is that there is a branch point singularity at v = +i which would 
imply a non-analyticity in the upper half plane, inconsistent with the definition of a 
retarded propagator. Note that this feature is purely a result of the high frequency 
approximation and the exact result (|3.35|) has no singularities in the upper half of the 
complex frequency plane. Indeed, closer inspection reveals that the putative branch 
cut originating at z> = +i has a vanishing discontinuity in the limit of large MR. The 
spurious branch cut originates from the equally spaced zeroes of Gr{v) appearing 
to coalesce in the high frequency limit. The branch cut discontinuity at v = — i 
is, however a genuine non-analyticity and originates from the infinite set of poles 
merging into a continuum in the high frequency approximation. The easiest way to 



Figure 6: For large mass and frequency, the propagator in the topological black hole phase, has a 
branch cut as indicated in the rescaled frequency plane, v — v/M. This results from the apparent 
merging of the infinite set of isolated poles of the exact Green's function. 
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understand the singularities and discontinuities of the function above, is to examine 
the function z 2MRz (z — i^- 2MR ( z - 1 ) an d then make the replacement z — > (1 — iP)/2. 

Along the branch cut B = {u G [—i, —ioo)}, we find that the retarded propagator 
has a discontinuity 



Disc| B G R (Mi>) =2i C A sia(2nMR) 



2 

z> = — z|z/| , \u\ > 1. 
It is clear that the discontinuity is large in the large mass limit. 

Now let us look closely at the putative branch cut along the imaginary axis 
B' = {D G [—*,*]}• Computing the discontinuity across this, we have 

'1 +x \ MR(l+x) /j_ xMJiM 



Disc | B ' Gr(Mu) = 2i C A sin (MRn(l - x)) 
v = ix, — 1 < x < 1. 



2 



(3.41) 

This vanishes when MR oo, for two reasons. The rapid sinusoidal oscillations 
will give vanishing contribution to any contour integral along B'. Furthermore the 
amplitude of the oscillations vanishes exponentially as evident from the expression 
above. 

In the leading high frequency approximation , for real frequencies, the imaginary 
and real parts of the Green's function are given by 

G R (u) « C A [^^) e-^T H tan m e ~2T H (cos(^) - i sgn(u) sin(^)) , 

T H = — , fel. 

(3.42) 

The high frequency, large mass limit thus appears to retain features of the thermal 
effects of de Sitter space. This result can be deduced from ( |3.40| ) after choosing an 
appropriate branch of the function and also directly follows from ( |3.37| ). We will see 
subsequently that these high frequency expressions can also be derived by solving 
the wave equations using a WKB approximation, providing a consistency check. 



3.5 R-current correlation functions 



Real time response to perturbations of a conserved charge density can reveal inter- 
esting late time physics, such as hydrodynamic or diffusive relaxation. It is now well 
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understood and established |28|, via holographic calculations in AdS black hole 
backgrounds, that correlators of conserved global currents exhibit hydrodynamic and 
diffusion poles. R-charge diffusion in the strongly coupled, high temperature J\f = 4 
plasma was first discovered in [P8| . The universal features of strongly coupled plas- 
mas follow from the properties of the stretched horizon of AdS black holes p9| . It 



is therefore natural to ask whether the horizon in the topological AdS black hole 
geometry implies hydrodynamic behaviour of correlation functions in the dual field 
theory. The answer to this question will depend on the relevant time scales involved 
since the boundary field theory is formulated on an expanding background, namely 
dS 3 x S 1 . 

The strong coupling correlators for 5*0(6) R-currents of M = 4 SUSY Yang-Mills 
theory on dS 3 x S l will be obtained holographically from the on-shell action for the 
SO(6)r gauge fields in the topological AdS 5 black hole background, by following the 
prescription of P3[. The Maxwell action for the gauge field is 



S = J d 5 x^-gg^g^F^F aPl (3.43) 

varying which gives the following equations of motion 

~^=f v {y/=ggr^Fap) = 0. (3.44) 

Here, we have gg G = I6n 2 R^s / N 2 . We find it convenient to use the following form 
for the metric in the region exterior to the horizon, 

ds 2 = Aids dz 2 + flL S (l-*) l dT 2 + CQsh 2 T dn 21 + T l d 2 ( g 45) 

Az 2 (1 - z) z 1 21 z 

Substituting the solution to the equation of motion that corresponds to the boundary 
value A a (2) „ = A° a back into the gauge field action, we will get a generating 
functional for the R-charge correlators of the field theory. 

Since the field theory lives on the dS 3 x S 1 boundary with spatial S 2 x S 1 spatial 
slices, it is natural to consider the late time behaviour of long- wavelength fluctuations 
in the following two cases: 



1. The R-charge perturbation is inhomogeneous on the S 1 , but homogeneous on 
the spatial section of dS 3 , 

2. The fluctuation is homogeneous on the circle, but inhomogeneous on the S 2 . 



Of particular interest is the presence or the lack of late time hydrodynamic relaxation 
of this system. 
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3.5.1 Inhomogeneous perturbation on the S 1 



In the first case, we will assume, for simplicity that the R-charge perturbation carries 
no momentum around the S 2 . Furthermore, we use gauge freedom to set the radial 
component of the gauge potential A z = 0. Hence, from the symmetries of the 
configuration, there are two remaining bulk gauge fields that are non zero: 

A T = A t (z,t,x), A x = A x (z,t,x). (3.46) 

Since the x-direction is a spatial circle, the two components of the vector potential, 
A T and A x , can be conveniently expanded in Fourier modes on the circle. The 
time dependence can also be re-expressed in terms of a mode expansion. There is 
a subtlety involved in this however. For A x , which is a scalar in dS$, the mode 
decomposition is straightforward 

in* poo J 

Axi*, ^x) = J2^J^T x (u, r) Q n (u, z), (3.47) 

where n G Z. On the other hand A T , which is a gauge field in dS^, has its complete 
time dependence captured by normal modes of two kinds - normalizable and delta- 
function normalizable. In fact, below we show that there is a single normalizable 
mode and a continuum of delta-function normalizable states obtained as solutions to 
a Schrodinger problem. Anticipating this we write 

p inx / r°° A,. \ 

A T (z, r, X ) = ^ [J ^ ^ ?M r) F n {v, z) + T r N (r) ^(z)j , (3.48) 

where T t (v,t) and T^{r) will be the delta-normalizable and normalizable modes 
respectively. The mode functions T X T are solutions to 



^cosh" 2 r-^- (cosh 2 rT r (z/,r))^ = -{y 2 + 1) T T [y, r), (3.49a) 



d 

dr 

cosh -2 r— (cosh 2 r T T ) = —i{y + i)T x 



The first of these can be put in the form of a Schrodinger equation, by defining 
T T = T / coshr, 

~f 2 -^f = u 2 f. (3.50) 

dr 2 cosh 2 r V ; 

It is solved by the associated Legendre function P 1 " w (tanhr). For u 2 > 0, there is a 
continuous infinity of delta- function normalizable states, which yield 

rr .. . e~ ivT /z/-itanhr\ ^ , . e~ ivT 2 
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The Schrodinger potential above also has bound states for v 2 < 0. In fact there is 
precisely one normalizable bound state with v 2 = —1, corresponding to the solution 

t " m ^Jt? (3 ' 52) 

This solution can be directly obtained by evaluating Pf^tanhr) at v 2 = —1 (and 
appropriately normalized) or can be systematically inferred from the Maxwell equa- 
tions. The continuum modes for A x and A T are orthonormal with respect to the 
inner product 

/oo 
dr cosh 2 TT a (is,T)T a (i/,r) = 2n6{is + i/), (3.53) 
-oo 

for a G {x, t}, while the bound state is normalized so that 

/oo 
dr cosh 2 r (T r N ) 2 = 1. (3.54) 
-oo 



In the far future the continuum modes are both given by T a ~ e~ lVT e~ 2r . Upon 
analytically continuing to the complex v plane, for v = i + u>, with w G 1, they 
are propagating (purely oscillatory) excitations with frequency uo in the far future. 
On the other hand the (normalizable) bound state decays exponentially in the far 
past and future and being real does not contribute to the flux at the horizon of the 
topological black hole. 



Using these modes to eliminate the r-dependence we find three equations that 
depend only on J- n {v, z) and Q n {y, z) 



i 



n AdS 

1 d 

^Ads dz 



, \ nil — z) 



o, 



n 



n 



Tn = 0, 



d , . 

Az Tz ^ 



zfG'n) -n(v + i)F n + (v 2 + l)G n 



0. 



(3.55a) 
(3.55b) 
(3.55c) 



Here, prime denotes a derivative with respect to z. The radial or z-dependence of the 
bound state solution, ^(z) is found by analytically continuing the profile for generic 
v to v — rki. Note that there are three equations for two unknowns; so, to ensure a 
non-trivial solution any two equations must imply the third and it is straightforward 
to check that this is indeed the case. We can then use these equations of motion to 
derive two independent ones, each containing only one of the unknown functions: 

Z/ 2 + l" 



4*(i - z)T': - 4(3z - - AT' n 
Az(l - z)QZ - A(5z 



n 2 - 



1-z 



8(1-2*) 

i)y n : — —y n 



n 2 



J ni 

(^ + 1) 



(3.56a) 
(3.56b) 



- 22 - 



where we have defined n = ni?AdS . These equations are immediately solved in terms of 
hypergeometric functions. Singling out the solutions that satisfy the purely infalling 
wave boundary condition at the horizon, we find the induced boundary action for the 
Maxwell fields (see Appendix A for details). The R-current correlators can be read 
off from the finite and non- analytic pieces of this boundary action ( |A.8| ), ([A. 12 ). We 



will denote the Fourier harmonics of the R-current along the spatial circle as 



2~ 



#(t)= / d X e- in xr(x,r), (3.57) 
Jo 

where we have already restricted attention to the s-wave sector on the spatial two- 
sphere. For perturbations with non-vanishing momentum along the spatial circle, we 
define the retarded, real time, Green's functions as 

CT{t, t'; n) = ([ j£(r), f_ n (r')\) 0(r - t>) , G { X , r}. (3.58) 

The conserved global currents are in one-to-one correspondence with the boundary 
values of the 5-d gauge fields. Functionally differentiating the induced boundary 
action (|A.12|) with respect to the boundary values of the gauge fields, we find 

G TT (T,T';n) = cosh 2 r cosh 2 rx (3.59a) 



n 2 



^T t (u,t)T t (-u,t') S(v,n) + T T N (T)T?(T') S(i,n) 

J-oo 271 



G xx (t, t'; n) = cosh 2 r cosh 2 r' [ — TJu, t)TJ-v, t') (is 2 + l) S(u, n), (3.59b) 

J-oo 2tt 

G TX {T,r';n) =G xr *(r,r» = (3.59c) 

cosh 2 r cosh 2 t n [ — t)T x (— u, t') (u — i) n), 
J-oo 2vr 



s (". n ) =Sf (I + 3(" ~v))+1>{\- s(" + ")))• ( 3 - 59d ) 

It is easily established that the above expressions are indeed retarded Green's 
functions and are non- vanishing only when r > r'. Two essential features ensure 
that this is the case: The function S(u, n) appearing universally in the //-integrals 
has only simple poles in the lower half complex plane at 

v = -i(2k+l)±n, keZ. (3.60) 

There is a second source of non-analyticities in the z/-plane. This lies in the in- 
dependent normalization of the mode functions T T (u, r) ( |3.51|) . The potentially wor- 
risome aspect of this is the appearance of a pole in the upper half plane at v — +i, 
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which gives a no n- vanishing contribution for r < r'. However, the potential prob- 
lem is eliminated by the term dependent on the discrete, normalizable mode T T N in 
( p.59a| ) which exactly cancels against the contributions from the poles at v — +i, 



ensuring that our Green's function is causal. The remaining Green's functions are 
manifestly free of singularities in the upper half plane. 

We thus see that the Son-Starinets recipe for determining real time response 
functions works in the case of the topological black hole, provided we carefully ac- 
count for the contributions from both the continuum and discrete mode functions in 
dS 3 . 

It is also clear in the above expressions, that there are no diffusion poles. Instead, 
from the properties of the digamma function which we have encountered before, the 
frequency space Green's function which is effectively S(f, n) , has only simple poles 
in the complex i/-plane, the lowest of these being at 

v — — i ± n. (3.61) 

Excitations with complex v = u — i, and u 6 R will propagate at late times, the 
perturbations simply evolving as left- and right-moving excitations on the S 1 without 
dissipating. 

The absence of any diffusion or transport like behaviour may be intuitively ex- 
plained by noticing the similarity of the topological AdS black hole to the BTZ black 
hole Setting up an excitation with momentum only on the S 1 is equivalent to saying 
the variation of the fields along the S 2 is zero. Therefore, aside from the time depen- 
dent factors associated to the de Sitter expansion, the metric that is "seen" by the 
bulk fields is not the full five dimensional metric, but its effective 2 + 1 dimensional 
portion, 

ds 2 = -R 2 AdS (p 2 - l)dr 2 + ^ [ dp 2 + p 2 r 2 x d X 2 (3.62) 

Comparing this with the metric for a 2 + 1 dimensional BTZ black hole with zero 
angular momentum 



ds 2 = - ( - Mj dt 2 + ^ - Mj dr 2 + r 2 dct) 2 , (3.63) 

we note the obvious similarity. Therefore we expect the behaviour of fields in the 
topological black hole background with an inhomogeneous excitation around the S 1 , 
to be similar to the behaviour of the fields in a BTZ black hole background. In other 



words, they should behave as in a (l+l)-dimensional CFT [27[|, just as we see from 
our results above. 
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3.5.2 Inhomogeneous perturbation on the S 2 



We will now examine the real time response to fluctuations carrying momentum 
along the spatial sections of three dimensional de Sitter space. Each spatial section 
of dSs is a two-sphere which undergoes exponential expansion at late times. For this 
case, we will focus on a situation where an inhomogeneous R-charge perturbation is 
set up on the two-sphere with only a dependence on the polar angle 9 and time r. 
For this configuration, the dual bulk gauge fields are 



A T — A T [z, t, 



A e = Ag{z,T,9); 



A Z = A X = 0, 



(3.64) 



where A z is set to zero by the gauge freedom, and A x vanishes due to ^-independence 
of the configuration. By spherical symmetry, the scalar potential A T and the vector 
potential Ag, each can be expanded in terms of scalar and vector spherical harmonics, 
respectively: 



A T = jVj(z,r) Y e °(9), A e = £>(*,r) dgY e °(9). 



Substituting these into the bulk Maxwell equations, we find 

4z(i - Z )d z ((l - z)d z g t ) - o 2 T g e + d T T, = o 



(3.65) 



(3.66a) 



4z(l - z)d z ((1 - z)d z Ti) 



^ {ft - d T Gl) = 0, 



cosh r 



(3.66b) 



d T (cosh 2 rd z Te) + i{i + l)d z G e = 0. (3.66c) 
^From equations ( |3.66b|) and (|3.66c|) we obtain a differential equation for T[ = d z Tt 



Ad z (z(l - z)d z ((1 - z)T' t )) - 



J~ p 



cosh T 



cosh r 



d 2 T (cosh 2 t\F;) = 0. (3.67) 



Now we will separate out the explicit temporal dependence, keeping in mind, as 
before, the possibility of contributions from both discrete and continuous modes 



/oo 
^r £ (z,,r)F^,^) + ^r £ N m (r)F; 
-oo , . 



(3.68) 



The mode functions % satisfy 



2 ^+!) 



-di- 



cosh r 



(cosh 2 r T(\ = v 2 (cosh 2 r TA 



(3.69) 
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which is a Schrodinger equation whose potential clearly will have both bound states 
and scattering or continuum states. The full set of solutions form an orthonormal, 
complete set. Indeed, the delta-normalized eigenstates are the Legendre functions 

rT .. , , . Pr^(tanhr) 

%[y, r) = T(l + iv) " \ 2 '-. (3.70) 

cosh T 

Those with v 2 > are scattering states with continuous values of v £ 1R, while the 
discrete, "bound states" have —iv = 1, 2, . . . £ 



_ m(i-m)\ Pr(tanhr) 
i<ra ~U (/ • m)! cosh 2 r ' ^ 7 j 



For i/ 2 > 0, the late time, r — > oo, behaviour of the modes will be significant, 



%(v, t) 



It>1 



e- WT e~ 2 \ 



(3.72) 



as these modes are oscillatory. Applying infalling boundary conditions on these 
modes at the horizon of the topological AdS black hole, ( |3.67| ) yields, 



F<,(v, z) = C t {v) (1 - z)- 1 -™ I 2 2 F X (-zf , 1 - if ; 1 - iv ; 1 - z) . 



(3.73) 



For the discrete series, the radial profile in the bulk, Ff m (z) is obtained by evaluating 
Fz(v, z) at v — —im. Putting the above ingredients together, the general form of the 
electric field along the radial direction in the bulk is 



1=0 

I 



dv . . P p ^(tanhr) 

— r(i + «/) -<— ^ — ^ f,(z/, * + 

Z7T cosh r 



m=l 



m (e + m V . cosh 2 r *<ml 



(3.74) 



This also allows us to automatically solve for A' e using ( [3.66c|) and we get 



A' e (z,T,i 



-E 

£=1 

e 



£(£ + 1) 



(l-m)! 



m 



(e+m)\ 



^ T(l + iv) <9 r P^(tanhr) F e (v, z) + 

Z7T 



a T Pr(tanhr)F £ N m (z) 



(3.75) 



m=l 



Now, the bulk gauge field action can be shown to induce a boundary term which 
will be the generating functional for the boundary R-current correlators. Using the 
explicit solutions above, the induced boundary action becomes 



S 



dr r\ 



J2 H*> T ) F' £ (z, r) cosh 2 r + £ £(£ + 1) Q t {z, r) g' t {z, r) 



1=0 



i=i 



z=e-*0 

(3.76) 
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The next step is to express this completely in terms of the boundary values of the 
gauge potentials .P°(r) = Tt{e,r) and G®{j) = Qi(e,T). Their radial derivatives 
T( and Q' e (equivalently A' T and A' e ) at the boundary z — e, are also determined 
completely by the boundary values of the gauge potentials, and G®{j) as in 



^From the boundary action above, we thus find that the real time, retarded 
Green's functions for the R-charge currents j^, in the gauge theory are 



G tt (t,t';£) =£(£ + 1) 



dv 7TV 



2tt sinh nv 



P < T^(tanhr)Pf (tanhr') T(i/)+ (3.77a) 



+ (-l) m m Pf(tanh r)P[ m (tanh r') T(z m) 



m=l 



G^(r,r';£) =£(£ + 1) 



C?Z/ 7TZ/ 



27T sinh7rz/ 



<9 r P - w (tanh r) <9 T , Pf (tanh r') T (i/) + 



(3.77b) 



+ (-l) m ma T P7 l (tanhr)o' T /Pf m (tanhT / ) T( 



2 m 



m=l 



G^(r,rV 



+ 1] 



27T sinli7rz/ 



P^* (tanhr) d T >P^ (tanhr) T(i/)+ 



(3.77c) 



+ ^(-l) m mP^(tanhr)^P^ m (tanhr / ) T( 



2 m 



m=l 



r(u) 



N 2 r 



W-f ) - &)) • 



(3.77d) 



We need to first confirm that these Green functions satisfy basic consistency checks. 
Specifically, the retarded functions must vanish for r < r'. As in the previous case, 
this property is not manifest, but follows from the nature of the non-analyticities 
of T(z/), and the normalized mode functions T(l + iu)Pf lu (tanhr), in the complex 



//-plane. For real values of u, the associated Legendre function is |4T|, 42 



r(l +iv)Pf iv (tanhr) = e~ iVT 2 F X (-£,£ + 1; 1 + iu, (1 - tanhr) /2). (3.78) 



For integer £, the hypergeometric function is a finite polynomial in tanh r and a ratio 
of degree £ polynomials of v. Therefore the exponential frequency dependence means 
that, for r — r' < 0, the integrals over the frequency v can be evaluated by closing 
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the contour in the upper half plane. The function T(z/) has no poles in the upper 
half complex plane. It has simple poles at 

u = 0, -2i, -4i... (3.79) 

The normalized modes, r(l+zi/)P £ ~ iI '(tanhr) have exactly I simple poles in the upper 
half plane at v = i, 2i . . . ii. The contributions from these are, however, cancelled by 
the inclusion of the discrete modes in the retarded Green's function above. Hence 
our correlators are zero for r < r'. 



3.5.3 Late time behaviour 



The real time response functions in general contain important information on the 
long time relaxation of perturbations away from the equilibrium or ground state. In 
thermal field theory on flat space, the relaxation of such fluctuations of conserved 
charges proceeds via hydrodynamic or diffusion modes. The response functions at 
strong coupling then exhibit diffusion poles in frequency space, G TT oc (iu — DA; 2 ) -1 , 
where u is the frequency and k, the soft spatial momentum. Due to the explicit 
time dependence of the background metric, we cannot do a similar frequency space 
study of the full Green's functions in de Sitter space. Instead, we could analyze their 
behaviour as functions of time. 



In de Sitter space, perturbations labelled by wave number £, get red-shifted so 
that given sufficient time their physical wavelengths become super-horizon sized. 
This happens when 

le~ T 1 

— ~ b; < 3 - 80 ' 

At late enough times, even very high harmonics on the sphere get stretched and 
eventually exit the horizon. To zoom in on the time evolution of such modes, it is 
useful to think of £e~ T , the physical wave number, as being fixed as r — ► oo. For 
example, in this late time approximation we neglect terms like I e~ 2r in comparison 
to powers of ie~ T . This is practically equivalent to going to planar coordinates for 
de Sitter space and the mode functions behave as 

P-(tanhr)|^ T=fixcd — r* J iv {2£e~ T ) (3.81) 

It is possible to derive this by replacing the potential £(£ + l)sech 2 r in the mode 
equation ( |3.69| ), with 4£ 2 e~ 2r . Note that, instead of a fixed physical wavelength 
if we focus attention on fixed comoving wavenumber, given by £, all modes simply 
approach the s-wave at late times, 

hm r(l + zz/)p--(tanhr)| £fixcd — > e*" . (3.82) 
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For fixed physical wavelengths, £e~ T , or equivalently, at the time when a har- 
monic I crosses the horizon, the real time correlators are given by the exact results 
with the replacement ( |3.81| ). The integral over v can be easily evaluated using the 
method of residues, and it turns out that the leading contribution at late times is 
from the residue at v = 0. Thus 



and 



G tt (t, 0; 



G 0t (t, 0; 



), G t °(t, 0; 



G*V,0; 



This late time behaviour is depicted in Figs. 
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Fig. 7: The leading behaviour of correlators G TT (and G T e), up to normalization constants: as a 
function of r on the left with I = 1000 and, on the right, as a function of le~ T 




Fig. 8: The leading behaviour of Gee (and Ge T )- The left figure is plotted as a function of time 
for I = 1000, while the right hand side figure is a function of le~ r . 

The late time behaviour deduced above is not characteristic of diffusion in de 
Sitter space. Suppose that the R-charge fluctuation relaxed via diffusion modes, then 
the covariant conservation of the R-current together with Fick's law would lead to 
the diffusion equation 

d T f = DV e V e ]\ (3.85) 

in dS%, D being the diffusion constant. The spherical harmonics of j T on the expand- 
ing spatial spherical sections would then obey, 

d T jI = -D £ -^ljl. (3.86) 
cosh r 
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At late times r — > oo and large enough £, this is solved by 

jJ~exp(|D^ 2 e- 2T ). (3.87) 



The large time behaviour of the Green's functions ( |3.83| ) and fl3.84|) do not match 



up with expected diffusive relaxation ( |3.87| ) on dS 3 . A natural reason for this is that 
the rate of exponential expansion of the spatial section and the Gibbons-Hawking 
temperature are both set by i? -1 . Thus the mean free path or the mean free time 
between collisions is comparable to the expansion time scales so that the system 
never enters a diffusive regime. 



4. The (small) AdS Bubble of Nothing 



The topological AdS black hole discussed above has a semiclassical instability when 

r, < f*f (41) 

which causes it to decay into a "bubble of nothing" in AdS space. The instability 
only occurs if fermions have antiperiodic boundary conditions in the ^-direction. 
With periodic boundary conditions for both bosons and fermions, the topological 
AdS black hole is absolutely stable. 

As originally pointed out in (and [|7| in the present context), the decay of a 
false vacuum in semiclassical gravity is computed by the Euclidean bounce which has 
the same asymptotics as the false vacuum in Euclidean signature. The bounce is a 
solution to the Euclidean equations of motion with a non-conformal negative mode. 
In the context of the asymptotically (locally) AdS spaces in question, the small 
Euclidean Schwarzschild solution represents such a bounce solution. In Lorentzian 
signature, the semiclassical picture of the decay process at time t = (say) involves 
replacing the t > part of the false vacuum solution (the topological black hole) 
with the appropriate analytic continuation of the Euclidean bounce to Lorentzian 
signature. The analytic continuation of the small Euclidean AdS black hole bounce 
which leads to dS 3 x S 1 boundary asymptotics, is the (small) AdS bubble-of-nothing. 
The metric for the AdS bubble-of-nothing solution is 



ds 2 = f(r) r\ d X 2 + fir)- 1 dr 2 + r 2 (-^- + cosh 2 (Jj-^ d^ 



f(r) = 1 + 



r 



1 rl(Rl dS + rl) 



(4.2) 



D2 2 
-"-AdS ' 
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Fig. 9: The global structure of the bubble geometry. The region inside the shaded region is empty, 
and its surface represents the de Sitter expansion of the bubble at r = r^ . 



In order to avoid a conical singularity in the interior, the periodicity of the compact 
X coordinate is related to as 



2nr x = 2irR 2 AdS " h . (4.3) 

' h "T -"-AdS 



Passing to the dimensionless coordinates introduced earlier, 

P = V(r/R Ad s) 2 + 1, r r h = ^, (4.4) 

the metric becomes 

r 2 _ n 2 

^2 _ jJ2 ,,,,, 'x , ff n \-l P J„2 i ( „2 t\r aJI i ™^2^r>2^ 



AdS 



f(p) d X 2 + f(p)- 1 ~^—r dp 2 + (p 2 - l)(-dr 2 + cosh 2 rd^ 2) 

K AdS P ~ 1 



(4.5a) 

Kp) =P 2 ~ ~TZ\ 'I {rl + 1) • (4.5b) 

In the AdS bubble of nothing spacetime, a slice of constant p is dSs x S 1 . The 
S 1 , however, shrinks to zero size smoothly at p = yjf 2 + 1 . The shrinking circle is 
the cigar of the Euclidean Schwarzschild solution. The boundary of the spacetime 
is approached as p — > oo. The semiclassical decay of the topological black hole at 
r = 0, results in the sudden appearance of a bubble of nothing in the region of 
spacetime, p 2 < 1 + r^. 
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4.1 WKB for the AdS Bubble of Nothing 



An exact holographic computation of correlation functions in the AdS bubble of 
nothing background appears difficult as analytical solutions to the wave equation in 
this background are not known. Despite this, we may obtain the boundary Green's 
function following a systematic approximation. In particular we will employ the 
WKB approximation to determine boundary correlation functions at high frequency 
and large mass. The WKB approximation has been used successfully |39| to find 



high frequency Green's functions in the Big AdS-Schwarzschild black hole. In this 
approximation, lines of isolated singularities (poles) get replaced by branch cuts since, 
in the high frequency regime, the separation between poles effectively goes to zero, 
as seen in our example above in Section 2.3.1. 

We first take the massive scalar wave equation 

d^V^g d u $) - M 2 $ = 0, (4.6) 



in the bubble background and expand the scalar field in harmonics on dS% x S 1 as 
in ( |3.3| ). The harmonics $ n (v, p) then satisfy a radial differential equation which can 
be put in the form 

d 2 $ n ( 1 1 1 \ d$ n 1 



dx 2 \x — 1 x + r\ x — r\ — \) dx 4(x — l)(x + rj^ix — r\ — 1) 



+ 7 + ^^ + 7 ~^ + T7 ,w . ~ow ~, 7TX 



(4.7) 



where x = p 2 and n labels the momentum along the S . This is an ordinary differen- 
tial equation with four regular or nonessential singular points at x = 1, — f^, r 2 + 1 
and oo. Analytical solutions for this type of equation are unknown. In fact, a sim- 
ilar differential equation was encountered in the computation of glueball masses at 
strong coupling in the three dimensional effective theory obtained from Euclidean 
thermal M = 4 SYM j5J (on M 3 x S 1 with SUSY-breaking boundary conditions). In 



that case the dual bulk geometry is the Euclidean black brane solution in AdS space 
where the thermal circle shrinks to zero size smoothly. 

The WKB solutions to the wave equation can be found after going to the 



- 32 - 



Schrodinger form by introducing the variables 

^Ttkry (4 ' 8a) 



u = — ^ cot- 1 (■?-)+ V1+r / coth- 1 I - F J== ) . (4.8b) 

These are the natural generalizations of ( ^.6| ) and ( |3.7| ) to the bubble of nothing 
geometry. The cigar in the geometry gets smoothly capped off at p = \Jr\ + 1, 
where the spacetime ends. In terms of the u coordinate, this occurs as u — > oo. 
In terms of the harmonics of \l/ on the dS 3 x S l slices, as in (|3.3| ), we obtain the 
Schrodinger equation obeyed by the harmonics ty n (y,u) in the (small) AdS bubble 
of nothing: 

du 




2 * n (i/, u) + V n (v, u) V n (v, u) = 0, (4.9a) 



p 2 -i;V p p 2 



K(^, «) = - -=— ^ _ 1 _ _ (rt + 1) + (4.9b) 



p 2 4p 2 



+ n 2 l^S + _^ (15p 4 _ 1Qp 2 _ i; 



Here p is implicitly a function of u, determined by the solution to (|4.8b ). In addition 



to the fact that the Schrodinger potential is far more complicated than ( |3.8bf) , one 
crucial difference to the TBH case is that the potential cannot be defined independent 
of the frequency itself. This is due to the term proportional to fh in (|4.9b|) . This 
situation is also in contrast to the case of the big AdS-Schwarzschild black hole in 
[39|. Here, for a given frequency (and mass), we need to find the "zero energy" 



eigenstate of the Schrodinger problem ( |4.9a| ). 

Interestingly, the qualitative behaviour of the potential V n changes, depending 
on the relative values of the mass and the frequency. This is illustrated in Figures 
(IUD and (ITTJ). The qualitative nature of the potential is easy to grasp in the high 



frequency limit, wherein we take 

v -> oo, MR -> oo, v= fixed. (4.10) 

MR y 1 

For simplicity, we also set n — 0, focussing attention on the spatially homogeneous 
fields. In this approximation, 

p ,4t) (> - 1 - 7 

The potential is vanishing (having thrown away a subleading constant in the large 
frequency limit) at the tip of the cigar p = ^Jr^ + l. Using u ~ — ln(p — y^f^-TT) 



V (u,u) - V(P,u) = (MR) 2 ( 1 - — — p 2 - 1 - - r\(rl + 1) ) . (4.11) 
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Figure 10: The Schrodinger potential for a massive scalar in the AdS bubble of nothing. In the 
above plot the dimensionless frequency v = 7, the mass MR = 2 and fh = rh/RAds = 1- 




Figure 11: The Schrodinger potential for v = 7, the mass MR = 7 and fh — ru/R^ds = 1- 



near this point, it follows that V(P, u) decays exponentially as a function of u. The 
potential has another zero at p = \Jv 2 + 1. Now, since the spacetime ends at p = 
y/l + rf l , the number of zeroes of the potential depends on whether \9\ is greater 
than or less than fh- In particular, when > fh, the potential energy has two 
zeroes or turning points, and qualitatively resembles Fig. |TU|, while for \u\ < fh, it 
behaves as in Fig. (|TTD (the potential asymptotes to a constant different from zero 
in the figures; however this constant becomes negligible in the high frequency limit). 



Zero energy Schrodinger wave functions in the potential ( [4.11|) are marginally 



bound states for |z>| > f h , since the potential has two turning points. For \u\ < f h , 
there is only one turning point and the wave function should exhibit a qualitative 
change in its behaviour at \P\ — fh- This also strongly suggests that in the complex 
frequency plane, v = ±f^ should be singularities of boundary correlation functions. 
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The zero energy wave function for the Schrodinger equation can be obtained in 
the WKB approximation, where care needs to be taken in applying the matching 
conditions at each of the turning points in the potential. We treat the two cases 
\u\ > r h and \P\ < f h separately. 

WKB approximation for v 2 > r\ 

There are two distinct regions in the potential: (i) the "quantum tunnelling 
region", y/1 + v 2 < p < oo which we label as Region I, and (ii) the "propagating 
region" , a/1 + f\ < p < \fl + v 2 which we call Region II. 

In Region I, we write the WKB solutions as 

*wkb(^,m)= y^j \ A+exp^J VV du^j + A_ exp ^— \ZVdttjj. (4.12) 

where the classical turning point u c is defined by 



p(u c ) = Vz> 2 + 1. (4.13) 

These represent the growing and decaying modes in the near boundary region of the 
bulk geometry. This can be understood easily as follows. 



du 
dp 



P 2 

p— >OC) " 



' =^ V w {MR) 2 p 2 (4.14) 



which then immediately yields the near boundary WKB solution (|4.12| ) for \E'. This 
together with its relation ( f4.8b| ) to the massive bulk scalar field $ implies 

*wkbUoo ~ A + (• • • )P~ 2 ~ MR + A -(---) P' 2+MR ( 4 -!5) 

where the ellipses denote unspecified normalization constants. The two power laws 
appearing in this solution are precisely the normalizable and non-normalizable modes 
of the massive scalar field, in the limit of large mass, in an asymptotically (locally) 
AdS spacetime. Following the prescription for computing the retarded Green's func- 
tions, we normalize ^wkb (|4.15|) so that it approaches unity near the boundary. 



In the interior, however, for p < yl + v 2 the solutions enter Region II and 
become oscillatory. In Region II, we have 



|V|V4 



^wkb(^,m) = TTTTT/i ( B+ ex P ( * / V M dp) + exp ( — i / v|V| du 



u 



"c 



(4.16) 

The constants are uniquely determined by the WKB matching conditions at the 
classical turning points of the potential V(v,u) and the one normalization condition 
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on A_ near the boundary. The details of the WKB matching conditions are explained 
in Appendix B. The crucial result of the matching procedure is that the solution in 
the near boundary region (Region I) has a normalizable mode with strength 



A, 



1 



A_ tan 



\/\V(u, u)\du 



(4.17) 



The argument in the above expession can, as usual, be identified with the action 
of a zero energy classical particle trapped in the potential V(v, u) between the two 
turning points u c and u — > oo (the latter corresponding to p 
spacetime ends). 



f\ where 



The crucial difference between the bubble of nothing geometry and the topo- 
logical black hole is the boundary condition imposed on the solutions to the wave 
equation in the interior. To extract retarded correlators from the geometry with a 
horizon we impose an infalling condition on the plane wave solutions near the hori- 
zon. In the bubble geometry, however, since spacetime ends smoothly in the interior 
where the cigar caps off, there is no freedom in choosing the boundary condition 
at the tip of the cigar - we must require regularity (normalizability) of solutions in 
the interior. This means that the solution to the Schrodinger equation ( 4.9a|) must 
approach a constant exponentially as u — ► oo. 




Figure 12: The exact numerical solution to the Schrodinger problem for the AdS bubble of 
nothing. Here vR = 200, MR = 100 and ru/R = 1. The solution approaches a constant for u ^> 1, 
while the classically forbidden region is < u < u c 0.49. 



The Green's function at high frequency is completely determined ( up to contact 
terms) by the ratio A + /A_. After substituting the solution ( 4.15|) into the boundary 
action we find, 

G R {p) « lim MR exp ^2 J y/V(i>,u) du^j tan y/\V(v,u)\ du^j . (4.18) 
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The exponential prefactor in this expression is the WKB transmission coefficient into 
the classically forbidden Region I. The physically relevant contribution to the trans- 
mission coefficient is the constant, e-independent term in an expansion around the 
boundary e — > 0. The leading e dependence, which is an overall multiplicative con- 
stant proportional to e 2MR , can be absorbed into the normalization of the correlation 
function. 



The first observation we can make without actually evaluating (|4.18|) , is that it 
an infinite f 
occur whenever 



has an infinite set of poles as a function of v on the real axis for v 2 > r\. These 



S n (p,MR)= / duy/\V{u,u)\ = (n+|)7r, n e Z. (4.19) 

J oo 



That is, the semiclassical action in the propagating region is a half-integral multiple 
of 7r. This is of course the condition for the existence of a bound state wave function, 
and the poles in the Green's function reflect the appearance of these bound states at 
certain values of v on the real axis. Recall that the corresponding correlators in the 
topological black hole phase do not have any poles on the real axis. 

At first sight the poles on the real frequency axis might appear somewhat sur- 
prising. However, they have the following natural physical interpretation: the low 
energy physics of the gauge theory on the boundary of the AdS bubble of nothing 
is that of nonsupersymmetric Yang-Mills theory at large N (and strong 't Hooft 
coupling) on dS%. The antiperiodic boundary conditions on the spatial S 1 make all 
fermionic excitations massive and the broken supersymmetry leads to large radiative 
corrections to the scalar masses. In the strongly coupled theory, the dynamical scale 
of the three dimensional effective theory is expected to be set by r^ 1 , the scale of the 
compact S 1 direction. When the Gibbons-Hawking temperature Tjj = 1/(2ttR) in 
dS 3 is smaller than r" 1 , we expect the gauge theory to be in a confined phase where 
the degrees of freedom are gauge singlet glueballs. The appearance of the isolated 
poles in the high frequency correlators is consistent with this physical picture. 

The WKB integral in Region II, in the high frequency approximation to the 
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Green's function (|4.18|) can be expressed in terms of complete elliptic integrals as 
S n {P, MR) = \MR 



dx 



V 2 

x-1 



2+1 ^(x + f 2 )(x-f 2 -l) 



MR 



V 1+2f l 



~2 ^ 2 

i/ >r h . 



(4.20) 



^From the general characteristics of these elliptic functions and their singularities 
41], it can be checked that Snip, MR) has no singularities on the real axis for 
> f\. Potential logarithmic branch points at v 2 = f\ and at v — 0, cancel out 
between the individual terms above. In fact, for any fixed value of fh, it also follows 
that, for large u, the WKB integral increases linearly with v. 



v 



Stt oc \u\ 



\i>\ > 1. 



(4.21) 



Hence for \u\ ^> fh, the propagator ( |4.18|) has approximately equally spaced simple 
poles on the real axis, whenever Su 



(n+|)7r. 



Although there are no other sources of singularities from Su, the WKB trans- 
mission coefficient in Region I, which also enters the Green's function ( |4. 18| ) can have 
branch point singularities on the real axis, 



S T (i>, MR) = —\MR 



dx 



1 - 



x-l 



1+P2 ^^ x + rl)(x-f 2 -l) 



MR 



V 1+2f l 



2c- 



\MR I In . 



(F (esc- 1 yfcD-K (f )) - i(l + rl) n (a|f) 
+ in 



(4.22) 



This function is also free of any branch cuts at v = ±f^, as can be checked by 
directly evaluating the integral at this point. However, the logarithmic growth at 
large v implies a branch point at infinity. 

As an aside we mention that the high frequency limit in the topological black 
hole phase ( |3.40| ) can be rederived by formally setting fh = in the WKB integrals 
and Gr{v) ~ exp(2«Sj). 
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WKB for v 2 < r\ 



For low (real) frequencies v 2 < r^, the nature of the WKB potential changes (Fig. 
(|TTD). Bound states are no longer possible. The zero energy WKB solution is 



*wkb(^,«) = A + exp (^J VVdu^j + A.^exp (- J Vv du^j . (4.23) 

The relation between the two coefficients is determined by matching to the wave- 
function at large u, where the potential decays exponentially and the wavefunction 
is a modified Bessel function (see Appendix B). We find that 



A. 



iA_. 



(4.24) 



The aymptotics of this solution near the boundary at u = e ( |B.26| ) allows to compute 
the boundary action and the Green's function 



G R (i>) » limi MR exp(<S 7 ) =iMR exp 2 / y/V(i>,u)du 



(4.25) 



In terms of elliptic functions the explicit form for the WKB action is 



A 



- r h 

X X X X • 



• — X X X x> » 



Figure 13: The analytic structure of the boundary Green's function in the " bubble of nothing" 
phase in the WKB approximation v — > oo, MR — > oo with v = v /MR fixed. Approximately equally 
spaced simple poles on the real axis for v 2 > ff l7 are accompanied by branch points at v — and 
infinity. 



s, 



MR 



V 



\(F{ 



CSC 



-1 



2e" 



\MR In . / r 



b) 

+ 1-K , . 



(4.26) 
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and the only singularity of this expression is the logarithmic singularity at v — 
which complements the singularity at v = oo found above. 

The resulting analytic structure of the Green's function is summarized in Fig. 
(|I3D, which is to be contrasted with corresponding Green's function (at large mass 
and frequency) in the topological black hole phase in Fig. (|5]). As already dicussed 
earlier, the isolated poles on the real axis indicate glueball states and that the gauge 
theory is in the confined phase, wherein the radius of the spatial S 1 is much smaller 
than the dS 3 radius of curvature. This may be interpreted as a hadronized phase, the 
de Sitter temperature being too low for the degrees of freedom to be deconfined. In 
this context, we should point out that the high frequency WKB analysis, where the 
frequencies are much larger than the de Sitter cosmological constant, is basically a 
flat space limit and thus the singularities of the Green's functions may be interpreted 
in the standard way as in flat space. One feature of the propagator in the bubble- 
of-nothing phase whose origin is not entirely clear is the branch point singularity at 
v = 0. A similar branch point at v = — i was encountered in the high frequency limit 
in the topological black hole phase. The associated branch cut was a consequence of 
the apparent merger of the infinite set of quasinormal poles of the topological black 
hole, in the high frequency limit. We do not know of a similar interpretation of the 
branch cuts of (ggg ) and Fig. (TfJ). 



5. Summary and Discussion 

In this paper, we have studied real time correlators in strongly coupled Af = 4 SUSY 
Yang-Mills theory on a time-dependent background, namely dS% x S 1 . In particular, 
we have calculated the retarded scalar glueball correlators and the R-charge current 
correlators in the Z^-invariant phase. 

The retarded scalar glueball correlators have an infinite number of poles in the 
lower half of the complex frequency plane, which represent the topological black hole 
quasinormal frequencies. The imaginary parts of these correlators are associated to 
the Gibbons-Hawking temperature due to the cosmological horizon of dS^. These two 
facts suggest that the Zn symmetric phase of the boundary field theory corresponds 
to a deconfined plasma in the exponentially expanding universe. 

We also computed the retarded correlators for the spatial spherical harmonics 
of the conserved R-currents using the Son-Starinets approach. Here we encountered 
a subtle point wherein we had to include in our mode expansions, the effects of 
real, normalizable, discrete solutions to the de Sitter mode equations, in order to 
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obtained a retarded Green's function. The corresponding frequency space correlator, 
appropriately defined in de Sitter space, also has an infinite number of poles in the 
lower half of the complex frequency plane, but these do not appear to correspond 
to diffusive poles. The lack of hydrodynamic behavior of the system is presumably 
due to the fact that the expansion rate of dSs is of the same order as the Gibbons- 
Hawking temperature. Here, we did not calculate the correlators of the stress-energy 
tensor. However, using the same argument as above, we expect not to find any 
hydrodynamic poles there either. 



In this paper, we have also calculated the retarded correlators of scalar operators 
Oa, with conformal dimension A ^> 1, in both the Z^-invariant phase and Z^ broken 
bubble phase. Unlike the correlators of the Z^r symmetric phase, the correlators in 
the Zjv broken phase feature an infinite number of poles on the real frequency axis. 
These poles are naturally associated to bound glueball-like states, which suggests 
that this phase is a hadronized phase, where the de Sitter temperature is too low 
to deconfine the degrees of freedom (the Hubble parameter is low compared to the 
dynamical scale of the effective, non-SUSY 3d theory on dS 3 ). Since this geometry 
contains no horizon, Son-Starinets prescription ^7f is not applicable in this phase 
and we have restricted ourselves to the high frequency and large mass regime by 
using WKB approximation. Since the prescription proposed by Skenderis and van 
Rees p5| , p6| does not rely on the existence of horizons in the geometry, it would be 
interesting to see whether one can obtain the retarded correlators beyond the high 
frequency limit using this prescription. 



It is also interesting to note that the relevant boundary condition on the horizon 
as prescribed by Son and Starinets |27J implies that the boundary theory is in the 
Euclidean or Bunch-Davies vacuum. This is in agreement with Ref. ||, where it was 
argued that due to the fact that the a- vacua Wightman functions for the topological 
black hole develop singularities on the event horizon, the preferred vacuum for the 
boundary field theory is the Bunch-Davis vacuum. It would be interesting to see 
whether one can understand the issue of a-vacua ambiguity better by applying the 
Skenderis- van Rees prescription |25|, |26| to this set-up. 
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Appendix A: Boundary action for bulk Maxwell fields 



A.l Nonvanishing S momentum 



We provide here, the steps in the calculation leading to the boundary action for 
fluctuations which are homogeneous on the spatial slices of dSs but with momentum 
along the spatial S 1 . The equations of motion ( |3.56a|) , ([3.56b|) can be solved in terms 
of hypergeometric functions 

T' n = C T {u, n) (1 - Z )-H"-W 2Fl (l + i(n-v),l-l(n + v);l-ip;l-z). 

(A.l) 

Here C T is a frequency dependent constant, to be determined by the boundary values 
of the fields. Near the boundary of AdS space, the solution approaches 



<p I _ (-1 r(i-w) 



(A.2) 



r^+|(n-,)jr^-|(n+„ 
x (2 7B + lnz + V(| + i(n-z/)) +^(i-f(n + z/))). 

Imposing the boundary conditions at z — > 

lim-F^e) = J*(v), HmG n {u,e) = £», (A.3) 

e^O e— >0 



from ( |3.55b|) we find that 

C T (u,n) = - ^ J}* L (n^ n {v)-*^n{v-i)Ql{v)). (A.4) 



With the normalization fixed in terms of the boundary values of the relevant fields 
we have that 

Kir, <0 =1 - n{u - i)Q°) x 

(2 7B + lne + V(| + §(n-i/))+^(i-|(n + i/))). 



^Frorn ( |3.55a| ) we also obtain 

= (A.6) 

We can now plug these solutions into the boundary action to obtain the retarded 
R-current correlators. Following identical steps for the normalizable (in time) modes, 
we have 

^'(e) = \n*F™ (2 lE + In e + if, (l + fn) + $ (l - in)) . (A.7) 
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The induced boundary action for bulk Maxwell fields has the form 

47T 



X 



I z=e 



(A.8) 

where </§ G = 16tt 2 R\dS / N 2 . The complete real time retarded correlation functions 
for the R-currents, can now be accessed readily. First we define the boundary values 
of our gauge fields as 

pinx 

A X:T (6,t, X ) = J2^ A xAt), (A.9) 



n 

so that 



/oo 
A n x (r)T x (u,r) cosh 2 r. (A.IO) 
-oo 

and similarly 

/OO fOC 
A n T (r) T T (u, r) cosh 2 r , ^ n N0 = / A n T (r) T r N (r) cosh 2 r . 
-oo </— oo 

(A.ll) 

Putting these ingredients together, we find that the boundary action is 



N 2 

S\ = —-— 

n 



/oo noc 
dr cosh 2 r / dr' cosh 2 r' 
-OO J —CO 

x (4"(r)<(r')fr x T T (i/,r)T T (-i/y) - ^Adsf [y ~ i) x 

(A.12) 

xT T (^r)T,(-,,r') - A-"(r)^(r') i? AdS f (1/ + i) T x (i/, t)T t (-^ r') + 
+^W^^f> 2 + l)r x (i/,r)T x (-i/y)) +r x A- T n (r)A n T (r')x 



x 



i£T n (t)^V) ( 27£ + ln e + ^ ^ + |^ + ^ ^ _ 



A. 2 Nonzero momentum along the spatial slices of dS 3 

Below we fill in the steps in the derivation of the boundary action for the Maxwell 
fields in the bulk. The asymptotic form of the radial dependence of the bulk potential 
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J 7 i, can be determined from ( |3.73| ) 

F / | z _ >0 = C / (i/)(C 1 (i/) + C 2 (z/) hxz) 



where 



Ci 
C 2 



2 7g + ^(i-f) + ^(-f; 

r(i-f)r(-f) 



r(i-f)r(-f; 



(A.13) 

(A. 14a) 
(A. 14b) 



We can solve for Ci in terms of the boundary values of the gauge potentials, using 
the bulk equation of motion ( 3.66b ) for T[ near the boundary which yields 

r°° flu 

4 / — r(l + ^)p-^(tanhr) G t {y) C 2 {v) = £(£ + 1) (J${t) - d T g° e (r)) . (A.15) 

Note that the other equation of motion ( |3.66a| ) for Q' e , just yields the time derivative 
of this condition, so that we have only one equation to determine the coefficient CV 
This equation can be solved if we recall that the associated Legendre functions are 
mutually orthogonal 3 



dr P; v (tanhr)P £ -^'(tanhr) 



5{v - v') 



T 1 -iv r (1 + iv) 



(A.16) 



Thus, we have 



Co 



4C 2 



-Y(\-iv 



/oo 
dr' (tf(r') - d T ,g e (r')) If (tanhr') (A.17) 
-oo 



from which we obtain the solution 

cosh 2 r^(e,r)= 
00 dv nv 



dr' {^(r')-d T ,g°(r')) 



2ix sinli7rz/ 



p-^(tanh r)Pf (tanhr') (lne + 2 lE + ^(-f ) + V(l - f )) 



m W$. f™(tanh T)if*(tanh /) (hi £ + 2 7 e + * (f ) + </>(! + f )) 



m=l 



(A.18) 



Plugging these back into the expression for the boundary action (|3.76|) , we obtain 
the generating functional for two point correlators of R-currents. 



3 The orthogonality of these functions for purely imaginary order follows from the 

fact that they are eigenfunctions of the Schrodinger equation in the sech 2 potential, 
d 2 



£(£ +1)/ cosh rj P\ v (tanhr) = v P™ (tanhr). In particular for v G R, these are scatter- 
ing states and are delta-function normalizable, and the eigenfunctions corresponding to two different 
eigenvalues are orthogonal as usual. 
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Appendix B: WKB matching conditions 



We explain below the matching conditions at the turning point (s) of the WKB po- 
tential for the Schrodinger equation (|4.9a ) in the AdS bubble of nothing background. 



B.l WKB matching conditions for v > 

Near the turning point u = u c corresponding to p = a/1 + u 2 , since we are well away 
from any extrema, we can assume that 

V(u, u) = k(u c — u) + . . . , u —* u c . (B-l) 

In this region where the potential is basically linear, the exact solution in terms of 
Airy functions is 

*U e = A + ^ Ai (k^K - u)) + A_ Bi («Vs (ttc _ u)) . (b.2) 

The normalizations and constants of integration have been chosen carefully so that 
the exact solution near the turning point, in terms of Airy functions, matches the 
WKB solution ( [4.12 ) in Region I (u < u c ) away from the turning point. Now, we 



can continue the solution ( p.2| ) into Region II (u > u c ). For u > u c , where the WKB 



solution (|4.16|) should be valid, the Airy functions have the asymptotic form 



*| sin(|v^( M - Mc )3/2 + i) eos(l^-^/ 2 + f) 

*\u>u c ~ Z/i + kV4( U - Uc )1/4 K l/4( u _ Uc )l/4 • 



Comparison with ( |4.16 ) near the turning point then implies 



5+ = e'i - ; 5_ = e^I (§A_ + iA + ) . (B.4) 

There is yet another condition that emergies from the behaviour of the solutions near 
the second "turning point" , p — > yl + f\ or w — > oo where the space ends. In this 
region we have 



I ^ a/ 1 + ^ 

d^Uy^Hf W - 2(2f2 + i)( p _ ^/lT^) 



(B.5) 



so that 



P~ « 2 V^ ex p(- 2 ^ M+ v% cot "(4?))- (B.6) 

It follows then that, as a function of it, the high frequency potential decays exponen- 
tially, 



VML__ » (Mi?) 2 ( 1 - ~2 J exp ( -2^=^ « + constants^ . (B.7) 
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Note that for v > r^, the potential approaches zero from below. Let us define 
constants A and B, in terms of which the potential is simply 

%teu)L M «-Ar*», (B.8) 

where A and B can be read off easily from the expressions above. The Schrodinger 
equation with an exponentially decaying potential is solved exactly by Bessel func- 
tions: 

- #"(«) - Ae- Bu V{u) = 0, (B.9a) 

= Ci J (2# e" Bu/2 ) + C 2 Y e~ Bu/2 ^) . (B.9b) 

Recall that we are looking for a zero energy eigenfunction of the Schrodinger prob- 
lem. This means that for a potential that vanishes at infinity, the corresponding 
(normalizable) wavefunction can only be zero. This is an important difference to the 
black hole case where the wave functions are infalling plane waves at the horizon. 
Requiring that the wave function \& vanish or approach a constant as u —>■ oo then 
eliminates the term proportional to Yq. Hence, in the exponentially decaying region 

oc Jo (2^e- Bu/2 ) . (B.10) 

The WKB approximation should match onto the Bessel function for large values of 
the argument of the Bessel function. Using the standard asymptotic expansion for 
Bessel functions 

j„WL,^^M + !^M. (B.ii) 



^From this we can deduce a relationship between the constants B + and B_ in ( 4.16j) . 
To make this precise we define the integral 



,r-\ / x 2 _ v 2 



S H (u) = I y/\V&Jf\du. = -MR I dxj {x2 _^ ){x2 + 1 + ^ (B.12) 



where 



S n (u(p)) = i 



MR 



i(l + fl) (n (a; sin"' JI£g\ kj - ^ H 



!^ 2 +l+f? , l+2f? , a 

a = o—^ ; b = -i h ; k = j. 

v- 1 ft ' b 



(B.13) 
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Then the WKB solution in Region II is 



WKB 



I V 1 1/4 



(B.14) 



For large u, 



J oo 



(B.15) 



Using this result and comparing (|B.14| ) to the asymptotics of the Bessel function 
( PUD , we find 

The final ingredient consists in determining A + and v4_. To this end we first define 



Si(u) = / \/V{y~u) du, 



(B.17) 



which then gives us 



Sr(u(p)) 



MR 



vF ( sin- 1 ( J\£$ ), k ) + 1 -(l + rl)U( a; sin 



,-i 



1 P 2 +rl 
a p 2 — l 



,1 ) 



(B.18) 



Near the boundary u — > or equivalently p — > oo, we find 



WKB 



^ 1^ p MR ~2 e inMR/2 



/MR 



l+r 2 h 



MR/4 



-MR/2 . 



( ' xp v^f W F ( cscr1 ^' fcJ-^-ja+^nW) 

/ A \ -MR/4 
_|_ ^ + 1^ p- MR -2 e -™MR/2 1 1 



exp 



'Mil 



1 + r? 



ly MR/2 x 



V 1+2f h 

Combining ( p.4| ) and (|B.16|) we obtain 



z/(F (csc^ 1 Vo, fc) - #(Jfe)) - ^(1 + f£) n(a|Jfc)J 



(B.19) 



= --4-2 tan(5//(u c )) . 



(B.20) 
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B.2 WKB matching for \v\ < 



When \v\ < r^, the potential energy V(v, u) is a monotonic function of u which 
exponentially vanishes as u — > oo. Now, the potential has effectively only one turning 
point and the wave function has no region where it propagates. The WKB solution 
( 4.12 ) in Region I should smoothly match onto the exact solution of 

-ip"(u) + Ae- Bu ip(u) = , A,B>0. 



(B.21) 



The solutions to these are the modified Bessel's functions. Enforcing regular be- 
haviour as u — > oo picks out 



The WKB approximation for Iq(x) is valid when i> 1, 

1 



lax) 



(e x + ie~ x ) . 



(B.22) 



(B.23) 



lx>>1 2V2™ 

We write the WKB solution to the wave equation in the bubble of nothing background 

as 



*wkb {y, u) = A + exp ( J Vvdu^j + A_ exp f- J Wdu^j 

Comparison with the modified Bessel function implies 

A + = iA— 

Near the boundary u — > 0, which is equivalent to p — > 00, we have 



(B.24) 



(B.25) 



WKB 



£ 1^ p MR ~2 e ™MR/2 



'MR 



1 + rl 



MR/A 



exp 



MR 



_|_ ^ + _L_ p- MR ~2 e -inMR/2 



exp 



1 MR 
MR 



,~2 



-MR/ A 



(B.26) 



z/ Mi ?/2>< 



>{F (esc- ^ fc) - ±K (i)) - 1(1 + f») ^B (a|| 
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